Ground state entanglement and geometric entropy in the Kitaev model 
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We study the entanglement properties of the ground state in Kitaev's model. This is a two-dimensional spin 
system with a torus topology and nontrivial four-body interactions between its spins. For a generic partition 
(j4, B) of the lattice we calculate analytically the von Neumann entropy of the reduced density matrix pA in 
the ground state. We prove that the geometric entropy associated with a region A is linear in the length of its 
boundary. Moreover, we argue that entanglement can probe the topology of the system and reveal topological 
order Finally, no partition has zero entanglement and we find the partition that maximizes the entanglement in 
the given ground state. 
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Entanglement is a topic of common interest in the two fields 
of quantum information and condensed matter theory. On one 
hand entanglement is the essential resource for quantum in- 
formation processing and quantum communication. Protocols 
based on entangled states can offer an exponential speedup 
with respect to classical computation Q]] . On the other hand, 
since entanglement measures quantum correlations for pure 
states, it is natural to think of a close connection between 
entanglement and the coiTelation functions of highly coiTe- 
lated states in condensed matter systems |2]. Entanglement 
seems to play an important role also in quantum phase transi- 
tions, where is believed to be responsible for the appearance 
of long-range correlations O01- Hence the quantification of 
entanglement in quantum states is a crucial problem for both 
fields. 

In quantum physics, it is important to determine whether a 
system can be simulated classically or not. Indeed, a quantum 
system that can offer a true quantum computation cannot be 
efficiently simulated by a classical dynamics. Recent works 
have demonstrated that any slightly entangled quantum sys- 
tem, i.e., a system in which the entanglement scales less than 
logarithmically with the number of qubits, can be efficiently 
simulated by a classical computer |6]. It is henceforth impor- 
tant to determine how the entanglement scales with the size 
of the system. To this end, we can define the geometric en- 
tropy \f\ associated with a pure state and a geometrical region 
A as the von Neumann entropy of a reduced density matrix, 
which is the measure of entanglement between the degrees 
of freedom inside the region A and outside it. The scaling 
of entanglement with the size of the system is connected to 
the holographic principle |8], i.e., the hypothesis that the en- 
tropy of a region is proportional to the area of its boundary. 
In the case of a critical spin chain, instead, the entanglement 
between a spin block of size L and the rest of the chain scales 
like S ^ log2 L and thus this system can be simulated classi- 
cally |4, 5]. 

In this article we study the ground state entanglement in Ki- 
taev's model [9J. This is a 2-dimensional exactly solvable spin 
system on a torus with four-body interactions in the Hamilto- 
nian. Its relevance stems from the fact that it was the first 
example of the new subject of topological quantum compu- 
tation 1I9I I10I1 and because it features topological order lllll . 
This is a type of quantum order that can describe those states 



of matter that are not associated to symmetries, like the frac- 
tional quantum Hall liquids. Topological order is not related 
to the symmetries of the Hamiltonian and it is robust against 
arbitrary local perturbations, even those that destroy all the 
symmetries of the Hamiltonian. 

By considering a subsystem of spins in a region A, we shall 
prove that the geometric entropy of A is linear in the length of 
its boundary. Moreover, we find that entanglement can probe 
the topology of the system and conjecture it can detect topo- 
logical order. 

The model. Consider akxk square lattice on the torus with 
spins-i attached to the links. There are 2fc^ links and the total 

Hilbert space H has dimension 2^*^ . The geometrical objects 
of the model are stars and plaquettes. A star is the set of four 
links sharing a common vertex. A plaquette is an elementary 
square (face) on the lattice. For any star s and plaquette p 
we define the star operator As and plaquette operator Bp as: 



As - n, 

system is: 
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The Hamiltonian of the 



H = -Y^As-Y^Bp 



(1) 



These operators share either or 2 links so they all com- 
mute: [j4s,i?p] = ,ys,p and hence the model is exactly 
solvable. Its ground state is £ — {\£,) G H : As\£,) = 
Bp\0 = 10' ^*iP}- We have the following two constraints 
on the stars and plaquettes: Yivs ^-s = IIvp ^p ~ ^' ^° '^°'- 
all the As , Bp are independent (there are only k^ ~ 1 inde- 
pendent operators of each type). Since they all commute, we 
can use them to label the states. By imposing the 2fc^ — 2 
independent constraints As — Bp — 1 we can label 2^'° "^ 
states in the ground state. Hence we find that the ground state 
is 4-degenerate, dim£ = 2^*^ "(^'^ ~^' = 4. More generally, 
the lattice does not need to be square. On a Riemann surface 
of genus g, the same model has a 4B-fold degenerate ground 
state, depending on the genus of the surface. This degeneracy 
is the sign of topological order 1 11]. 

Given a curve 7 running on the links of the dual lattice, we 
define the associated string operator iy^[7] — Yijey'^V-' ^^ 
j e 7 we mean all the links crossed by the curve 7 connecting 
the centers of the plaquettes. A string-net is a product of string 
operators. 



Let A be the group generated by the /c^ — 1 independent star 
operators Ag', \A\ = card^ = 2*^"^^. The elements g ^ A 
are products of closed contractible strings |9,il3]. Let 71, 72 
be the two non-contractible loops generating the fundamen- 
tal group of the torus tti{T'^) and define the associated string 
operators ^1.2 = ^^^[71,2]- We refer to Wi as "ladder" oper- 
ators since they flip all the spins along a ladder going around 
a meridian, respectively a parallel, of the torus. Let N be 
the group generated by A,'Wi,W2- This group contains all 
the possible products, contractible or not, of closed strings. 
Closed strings and closed string-nets commute with H, while 
open strings do not commute j^. The group N splits in four 
cosets with respect to ^ : Ae = Noq,Awi = Noi,Aw2 = 
Nio,Awi'W2 = Nil. Since N is Abelian, right and left 
cosets coincide. The four cosets are the equivalence classes in 
W/A and all have the same order 2'' ^^; hence \N\ = 2*" +^. 

Let S be an orthonormal basis of the Hilbert space H, S ^ 



{|e> 



,S2kV 



0, 1}, where each Sj labels the 2- 



dimensional Hilbert space of a single spin/link in the lattice. 

Let |00) = \si = 0, ...,S2/£2 = 0). Consider the orbit of 
|00) through iV, S' = N\00) = {g\00)\_g e N}. Thus S' is 
the homogeneous space with respect to N. Since the action of 
A^ on 5 is free, we have \S'\ — 2^ +^ 

The ground state. We now give an explicit construction 
of the vectors |^) in the ground state L. The constraint 
-SpIO ~ 10 7 ^P means that the^round state is a superposi- 



E 



|ej)e5' "JPJ/ 



tion of only basis vectors in S' |9], \^ 
Egelva(.9).9|00). 

The constraint As|0 = 10'^* implies that all the vec- 
tors from each sector iVy |00) are in the superposition with 
the same coefficient, that is |i^) € £ if and only if |^) = 
l-4r'/'EL=o«^^.E,6lv.,5|00) withEL=o|a,,f = \. 
This is because ^s E^elv,, 5|00) = Egelv., 5|00> since 
each sector is mapped into itself by star operators, since As E 
N, AgNij = Nij and there are no invariant subspaces within 
Nij. We can thus write the ground state subspace as £ = 
span{|C.,), I, J = 0,1}, where |6,) = \A\-^'^Y.geA9M 
and \ij) — w\w\\00). By construction |^y) are orthonormal. 

Entanglement properties of the ground state. Multipartite 
entanglement is notoriously difficult to quantify, since there 
is no known entanglement measure for the general state of a 
quantum system. However, for an arbitrary pure state pab 
of a bipartite system (A, B) there is an (essentially) unique 
entanglement measure, the von Neumann entropy S: 



S = -Ti{pA log2 pa) 



(2) 



where pA = TiBiPAB) is the reduced density matrix of the 
sub-system A. 

Proposition 1. For a given lattice partition {A, B), the four 
ground states |^oo), ICoi). |?io), iCii) have the same enti'opy 
of entanglement S. 

Proof. Let Wi = Wi^A ® Wi,B, where Wi^A(B) acts only on 
the A (B) subsystem. From the circular property of the trace 



and wf = 1, w 



1, 2, it follows immediately that all 



four ground states have isospectral reduced density matrices, 
e-g-, Pa{£,oi) = Trs(|^oi)(^oi|) = Trs(u;i|^oo)(^oo|«^i) = 



wi^aTtb ( 1^00) (Coo Dwi^A- Therefore S is the same for aU four 
ground states. D 

Any element g E N can be decomposed asg~XA^XB, 
with xa{xb) a product of spin flip operators J| a^ acting 
only on the subsystem A{B). In general xa.b are not neces- 
sarily closed string nets. If, however, xa is a closed string-net 
(from which follows immediately that also xb is), we will 
write g = gA® gB, gA,B e N. 

We introduce now the following subgroups of A acting triv- 
ially on the subsystem A and respectively B: 



Aa = {g^A\ g = gA® Is} 
Ab = {g(^A\ g^^A® gB) 



(3) 
(4) 



The order of these groups is dA = \Aa\ and ds = l^lsl. De- 
fine also the quotient A/ Ab and let / = |^/^b| = I^I/^b 
be its cardinality. Notice that / is the number of elements in 
A that act freely on A. If there are n independent stars As 
acting on A, it turns out that / = 2". In other words, / is the 
number of different configurations of spins on A compatible 
with the states in S' . Define the group Aab = A/ {A a ■ Ab)- 
Wehavel^AsI = \A\/dAdB- 

Lemma. The reduced density matrix corre- 

sponding to the basis vector |^oo) is Pa = 

/"^ Y.geA/AB XA\OOA){OOA\xAgA- 

geAA 
Proof. First, recall that any two group elements 

g,g' € A are related by g' = gg, with g € 

A. We can write the reduced density matrix for the 

state ICoo) as pA = 1^1"^ Tr^ ^^ ,,,g_4g|00)(00|<?' = 
\A\~^Y.g,g'eA^A\OOA)(OOA\x'A{OOB\x'BXB\OOB);lhenon- 
zero elements of this sum satisfy x'^xb = 1, which im- 
pHes x'^ = XAgA, with g = cja ® ^b € Aa- There- 
fore pa = l-4|"^Ege^,geyt^^'4|00A>(00A|a;A5A- Now, 
x'ji = Xa if and only if g' — hg, with h e Ab, and since 



/ = ds/l^l, we can write 

PA = f-' Y. ^a\00a){00a\ 

qGA/Ab 
S^Aa 



XAgA 



(5) 



Notice that if dA = 1, the reduced density matrix is 
diagonal. D 

We are now ready to prove the main result. 

Theorem. Consider a partition {A, B) of the lattice, and 
let the system be in the ground state |Coo)- The entropy of 
entanglement is S = log2(.//rfA) = log2 \A\—log2{dAdB) — 

log2\AAB\- 

Proof Let us compute the square of the reduced density 
matrix: 



P\ 



= /-2 J2 Xa\00a) {00 A\xAgAx'A\00 a) {00 A\x'A~g'A 



g,g'£A/AB 
g,g'eAA 



= /-2 J2 Xa\00a) {00 A\xA~gA~g'A 



geA/As 
g,g'eAA 



f-^\AA\ J2 XA\OOA){OOA\xAgA = f-'dAPA (6) 



geA/AE 
g&AA 




FIG. 1: A region A of the lattice obtained by taking all the spins 
(solid thin lines) inside or crossed by a convex loop (thick line). The 
length LeA of its boundary is equal to the number of spins intersected 
by the loop and this is equal to the number of (white) sites on which 
the stars act on both A and B. The area Ea of A is equal to the 
number of (gray diamonds) sites inside it and cIa ~ 2^-^ . Then the 



area of B (black sites) is Eb = k^ 
entropy is S = LaA — 1. 



Sa — LgA- The geometric 



From the Taylor expansion of the log we obtain logj pa 
PAJd'^ log2(dyi//) and the entropy of entanglement is 



S^\og^{fd-/) = \og^ 



\A\ 
dAd-B 



l0g2 1-4 



AB 



(7) 



The meaning of this formula is the following. Both Aa and 
Ab contain information related exclusively to the subsystems 
A and B, respectively. Hence we say that it is the informa- 
tion contained in the bulk. Then the entropy is given by the 
difference between the total disorder and the disorder in the 
bulk. This implies that the entropy is due to the boundary 
between the two subsystems. If we choose the partitions in 
a convenient way, we can give a clear geometrical picture of 
this boundary. Consider now A as the set of all spins inside 
or crossed by a loop in the dual lattice (see Fig[0. The spins 
intersected by the loop are the boundary of A, while the ones 
inside are the bulk. The area E^ of A (in lattice units) is the 
number of sites inside the loop and dA — 2^-*, since this 
is the number of independent closed string-nets acting exclu- 
sively on A. Similarly cIb — 2^^, where E^ is the number 
of star operators acting exclusively on B. Then the entropy is 

S = l0g2(2'='-l) - log^idAdB) ^P~l~^A~^B. If 

the loop is convex (i.e., a rectangle), the perimeter LgA of A 
is equal to the number of spins intersected by the loop and this 

is LgA = fc^ — J^A — ^B (see FigQ. Therefore 



S — LgA — 1 



(8) 



Notice that in this case the entropy of a generic ground state 
1^) is also S = LgA — 1. This follows from the fact that the 
ladder operators can always be moved outside the disk (see 
also til). 

If now we take as the boundary of A an arbitrary (con- 
tractible) loop 7 on the lattice, the entropy will be 5 = 
ai{j)LgA — 1 + 02(7), where ai(7)'s depend only on the 
geometry of the loop II2I1 . Hence in general the entropy is 



linear in the boundary length of A. We see that this entropy 
is the geometric entropy |7] obtained by computing the pure 
state density matrix and tracing over the variables inside the 
geometric region A to obtain a reduced density matrix to eval- 
uate the entropy S. We want to explain the very important fact 
of the presence of the —1 in Eq.(|S}. The —1 comes from the 
constraint that we have on the number of independent stars; 
this follows from the fact that each link is shared by two sites, 
hence Ylvs ^s = I- This construction is not restricted to a 
square lattice. On a Riemann surface of genus g the number 
of sites, links and plaquettes {no,ni and ^2, respectively) is 
constrained by the Euler's formula: uq — ni + n2 — 2(1 — g). 
The number of independent stars and plaquettes is again no — 1 
and n2 — l (the last constraint follows from the boundary con- 
ditions); hence we can label rig + 77,2 — 2 = ni — 2g states 
and thus the ground state is 2^B-fold degenerate. The system 
exhibits topological order lllll . Now we can sketch an argu- 
ment to compute the entropy for the surface of genus g. If the 
number of stars and plaquette is the same (no = n2), then the 
entropy is 5 = ni/2 — g — E^ — E^. The topological order 
in this model manifests itself in both the degeneracy and the 
entanglement in the ground state. This suggests the very ap- 
pealing possibility that entanglement could detect topological 
order 

In the following we calculate the entropy S for several sub- 
systems A of the lattice. From Proposition 1, 5 is the same 
for all the states |^y), hence we will consider the system to 
be in the |^oo) ground state. For a generic ground state |^) the 
results will be presented elsewhere \VZ\ . 

Example 0: one spin. If the subsystem A is a single spin, no 
closed string net can act just on it; hence dyi = 1, / = 2 and 
the entropy is S* = 1. This shows that any spin is maximally 
entangled with the rest of the system. 

Example L two spins. We now compute the entangle- 
ment between two arbitrary spins on the lattice. In this 
case we use as an entanglement measure the concurrence 
C. For a general mixed state pij of two qubits, the concur- 
rence is C = max{0, ^/\i — \/A2 — ^f^ — ^/\i\, where 
Ai, A2, A3, A4 are the eigenvalues (in decreasing order) of the 
matrix pijia^ ® ay)p*j{ay ® a^) ijlj. It is easy to see that 
no g G A acts exclusively on the two spins, so d^ = 1 and 
the reduced density matrix pij is diagonal (we considered the 
two spins as the A subsystem). Let pij = diag(a, b, c, d); in 
follows immediately that C — always, hence there is no 
two-qubit entanglement between any two spins. 

Example 2: the spin chain. We take the partition A con- 
sisting of all the spins belonging to a meridian (or parallel) 
71 of the torus; this is a system of k spins-i. It can be eas- 
ily seen that no closed string-net can act exclusively on A and 
thus d^ = 1. The number of possible configurations of spins 
on the chain 71 is 2^, but there are only / = 2^^^ configura- 
tions of spins in A that enter the ground state, namely the ones 
with an even number of spin flips. We see indeed that we have 
fc — 1 stars acting independently on the chain and we can ob- 
tain all the allowed configurations applying products of these 
stars (i.e., elements in A acting freely on the chain), which 
gives / = 2^^^. Then the entropy is 5 = log2 / = fc — 1. 

Example 3: the "vertical spins ". Let A be the set of all ver- 



tical spins on the lattice (and hence B is the set of all horizon- 
tal spins). This partition is such that no closed string operator 
can act trivially on either subsystem, hence (Ia = ds = I and 
we obtain the maximum possible entanglement S = k^ — I. 

It is interesting that no partition has zero entanglement, 
since dAds < \A\ for any (nontrivial) partition of the lattice. 

Conclusions. Entanglement is one of the most striking fea- 
tures of quantum mechanics and the most important resource 
for quantum information processing. Entanglement between 
two parts of a quantum system is measured by the von Neu- 
mann entropy. The entropy of entanglement can be connected, 
with some assumptions 1 14], to the notion of entropy in ther- 
modynamics. A whole number of studies is devoted to inves- 
tigate how the entropy associated with a region is related to 
the geometric properties of that region. A simple argument 
by Srednicki [151 for a free massless scalar field says that en- 
tropy should scale like the area of a spherical surface enclos- 
ing the field. This is because the entropy associated with the 
volume inside the sphere must be equal to the entropy asso- 
ciated with the volume outside the sphere. Then the entropy 
should depend only on properties that are shared by the two 
regions (inside and outside the sphere), that is, the area of the 
shared boundary. Srednicki |15] has calculated (numerically) 
the entropy for 2 + 1 and 3+1 conformal field theories. It 
turns out that the entropy associated with a region TZ is pro- 
portional to the size ^(T?.) of its boundary. This result is in 
agreement with black-hole thermodynamics, whose entropy 
is Sbh = MpiA/A, where Mpi is the Planck mass and A 
is the surface area of the horizon of the black hole flol- A 
recent work 1 17] proves analytically that the entropy associ- 
ated with a region of a discretized free quantum Klein-Gordon 
field can be bounded from above and below by quantities pro- 
portional to the surface area. Several authors have calculated 
the entanglement in ID spin chains. In the case of XY and 



Heisenberg models, the authors in Refs. f^ have calculated 
the entanglement between a spin block of size L and the rest 
of the chain. They found two characteristic behaviors. For 
critical spin chains, the entanglement scales like S ^ logj L, 
whereas for the noncritical case S saturates with the size L of 
the block. This result is in agreement with the result for black- 
hole thermodynamics in 1h-1 dimensions flo.M^, which sug- 
gested a connection between the entanglement measured in 
quantum information and the entropy of the vacuum in quan- 
tum field theories. 

In this article we have investigated the ground state entan- 
glement in the Kitaev's model. This is a two-dimensional spin 
system with four-body spin-spin interaction which presents 
topological order. On a Riemann surface of genus g the degen- 
eracy of the ground state is 4^ and such degeneracy is stable 
against local perturbations. This is at the root of the topolog- 
ical quantum computation. We found that although no two 
spins of the lattice are entangled (the concurrence is zero for 
any two spins), the ground state has multi-body entanglement. 
For a generic partition {A, B) of the lattice we calculated ana- 
lytically the von Neumann entropy of the reduced density ma- 
trix PA in the ground state |^y ). We found that the geometric 
entropy associated with a region A is linear in the length of its 
boundary. Moreover, no partition gives zero entanglement, so 
the system has an absolute entanglement. Finally, we argued 
that entanglement can probe the topology of the system and 
raised the very interesting question of whether it could detect 
(or measure) topological order 

RZ. gratefully acknowledges funding by European Union 
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Note. After the completion of this work, we came across 
Ref. 12011 which uses a similar formalism in the context of 
stabilizer states. 
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